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We briefly review some of the developments in the study of parton distributions 
which have occurred since DIS2000, including discussion of uncertainties, shadow- 
ing, unintegrated and generalized distributions. 



1 Data and parton distribution functions 

The situation is summarised in Fig. 1, which shows the kinematic regions 
in the (x, Q 2 ) plane covered by (i) the experiments at HERA, (ii) the fixed 
target deep-inelastic scattering experiments and (iii) the single jet inclusive 
experiments at the Tevatron. Roughly speaking, the fixed target experi- 
ments determine the distributions of the u,d,u,d,s,s quarks for x ~ 0.1 
(and Q 2 ~ 10 GcV 2 ). The HERA measurements determine the sea quark 
distributions, and their dF^jd In Q 2 data determine the gluon, for x ~ 10~ 3 
(and Q 2 ~ 10 GcV 2 ). The Tevatron jet data determine the quarks and gluon 
distributions in a region around x ~ 0.1 (and Q 2 ~ 10 4 GeV 2 ). 

In the past year new data have been available in each of the three domains. 
The HI a and ZEUS 13 experiments have measured Fj_over a larger domain 
(see Fig. 1) with much improved precision. The DO and CDF □ collabo- 
rations have measured the inclusive single jet Et distribution and provided 
correlated statistical and systematic errors. In particular, the DO experiment 
has measured the Et distribution in five different rj bins, out to rj = 3 and so 
samplesjpartons in a much wider range of x (see Fig. 1). The NuTeV collab- 
oration have used is, v beams to measure F2,xF$ with increased precision, 
and also, by observing fj, + production, to obtain information on the s and 
s quark distributions; CCFR have made a model-independent re-analysis El of 
their v, v data, which removes the discrepancy for x < 0.1 with NMC. The 
E866 collaboration a have observed both pp and u pn" Drell-Yan production 
and further constrained the difference of the u and d distributions. 

2 Uncertainties in parton distribution functions and 
observables 

The parton distributions are determined from NLO DGLAP fits to data from 
about 14 diverse experiments, with typically 1500 data points and more than 
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Figure 1. The kinematic domains probed by the various experiments, shown together with 
the partons that they constrain. The new data that have become available since DIS2000 
are also indicated. 



20 parameters. There are many sources of uncertainty. First the statisti- 
cal and systematic experimental errors. The latter are often not randomly 
distributed and, moreover, may depend on theory. The theoretical uncer- 
tainties come from higher order QCD contributions, the choice of factoriza- 
tion and renormalization scales, resummation corrections (from resumming 
lnl/x,ln(l — x) terms), power law contributions, nuclear target corrections, 
the choice of the particular parametrization of the starting distributions and, 
in hadronic processes, from the treatment of the underlying event. 

To date, the best attempt to determine the errors on parton distributions 
is due to Botje 0. He includes the statistical and systematic covariance matri- 
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ces and allows for some of the theory uncertainties. However, the analysis is 
based on five massless flavours and so is only applicable for x > 10~ 3 . More- 
over there is no constraint on the gluon at large x values (as would come from 
including Tevatron jet data in the fit) and so the gluon is undetermined for 
large x, which feeds down to small x via the momentum sum rule. A similar 
analysis, also based on the Hessian method, has been performed by CTEQB. 

Recently there have been several contributions concerning the errors on 
partons and observables. Giele et al. a have obtained 100,000 'optimized' 
sets of partons, and expressed the uncertainties as a density measure in pdf 
functional space. In principle it is 'easy' to propagate errors given complete 
data information, but the question of theory uncertainties has to be addressed. 
So far only data have been considered; Giele et al. conclude that NMC 
and (old) ZEUS data are incompatible with the other F^ 9 data sets. This 
type of analysis is in its infancy. Clearly it is necessary to include other types 
of data, e.g. neutrino and jet data. 

It is advantageous to determine the uncertainties on observables directly, 
rather than to go through the intermediate stage of determining the parton 
errors. The disadvantage is that global analyses have to be carried out for 
each observable separately. In view of the importance of W and Z boson 
hadroproduction as luminosity monitoring processes .at the Tevatron and the 
LHC, attention has focused on aw,z- MRST(1999) LL3 found an uncertainty 
Aaw of ±3% at the Tevatron and ±5% at the LHC. In the latter case the 
uncertainty is mainly due to the-/ conservative) error of ±0.005 assigned to 
as(Mf). More recently CTEQ O found Aa w = ±4% at the Tevatron and 
±8% at the LHC, taking only data errors into account. These surprisingly 
large errors are due to using older data and choosing a large tolerance, \ 2 ~ 
Xmin < 1^0. In a neaz-global analysis, which incorporates all the new precise 
data, MRST (2001) B find Aa w = ±2% at both Tevatron and LHC energies, 
when only data errors are considered. However theory errors (Aas etc.) have 
to be included. Interestingly this latest analysis, which incorporates the new 
Tevatron jet data EKa, gives an error of only ±20% on the high x gluon density. 

3 The description of F 2 data at low x 

A good description of the F% data in the low x domain can be obtained in 
many different ways. First, there axe empirical models, with very few input 
parameters (such as that of Haidt E3) , which give excellent descriptions down 
to very low Q 2 . Then there are dynamically-motivated parametrizations, such 
as the 'saturation' model of Golec-Biernat and Wiisthoff EJ based on the qq 
dipole framework. Surprisingly, on a more comprehensive level, the NLO 
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DGLAP analyses Ejl3'LH13^y^t3 continue to give satisfactory descriptions of the 
now very precise F% data, down to remarkably low Q 2 and x. However in these 
fits the gluon tends to go valence-like or negative at low Q 2 , Q 2 ~ 1 GeV 2 , 
indicating that the approach ceases to be valid in this domain. This is also 
reflected in the anomalous behaviour of Fl in this domain, perhaps indicating 
the need for ln(l/ir) resummations. In fact, a good description of F2 is also 
obtained in a unified approach which incorporates DGLAP and BFKL with 
higher order ln(l/x) contributions 113. Here a 'flat' gluon is input. 

The message is clear. One low x observable can be described in many 
ways. It is possible to trade effects of the different types of perturbative 
evolution with different choices of the input forms. Of course, the situation 
would be changed if another independent quantity, such as Fl , were precisely 
measured at low x. 



4 BFKL and DIS 

The resummation of the ln(l/ir) contributions is now known to NLO El. The 
resulting x — > behaviour has the form x~ u with 

u = u (l - 6.5as), (1) 

where w = <5s4 In 2 is the LO behaviour; as usual as = 3ag / tt. At first it was 
thought that such large NLO corrections would mean that no stable small x 
predictions could be made using the BFKL procedure. However it turns out 
that this is not the case. It is possible to identify higher-order terms and 
then to resum them. Indeed Ciafaloni et al. t2l carry out an all-order ln(l/a;) 
resummation of the following effects: (i) running as, (h) the non-singular 
.DGLAP terms and (iii) the angular ordering and energy constraints. (Thorne 
E3 has made a recent detailed study of effect (i)). The result is a stable x~^ 
behaviour which is consistent with observations. In fact, prior to this, the 
fit of Kwiecinski et al. E-3, mentioned at the end of Section 3, was based on 
a unified equation which incorporates these all-order ln(l/a;) contributions, 
where the imposition of a consistency (or kinematic) constraint Eil plays a 
major role. 

We may use a very simplified calculation to show how including effect 
(iii) tames the NLO behaviour of (]l|). Recall that as we proceed along the 
BFKL gluon chain, we have ordering in the longitudinal momenta, Y > Y' , 
where Y = ln(l/a;). To allow for the energy/angular ordering constraints, the 
ordering takes the modified form 

Y > Y' + S. (2) 
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Hence the BFKL equation for the unintegrated gluon distribution becomes 

fQO = fo + I Kf(Y')dY', (3) 
or, in differential form 

df(Y) 
dY 



wof(Y-S). (4) 



Note that if S — then / = f exp(uj Y) ~ x u ">, as required. In general, the 
solution of (Q) is / ~ exp(wF) ~ x~ u , where 

lu = uu e - """ 5 = ^o(l-^o^ + •■■)■ (5) 

Comparison with (|l|) reveals S = 6.5/4 In 2 = 2.3, so 

u) = luq exp(— 2.3a;), (6) 

which-Lcads to a similar behaviour of lu as a function of as to that found in 
Refs. This toy model illustrates how the inclusion of a summation of 

higher-order terms stabilizes the NLO result. 

A more 'phenomenological' way of performing the resummation of higher- 
order contributions, and achieiiing the perturbative stability of the BFKL 
approach, has been proposed E3. This incorporates the one- and two-loop 
terms in the anomalous dimension, imposes momentum conservation and 
parametrizes the residual ambiguity in terms of a single parameter LO Q g which 
specifies the x~" off behaviour as x — * 0. 

5 Parton shadowing 

Let us start with the original GLR equation El, written in terms of DGLAP 
evolution for the gluon density 



d(xg(x,Q 2 )) _ N c a s 



x 9 N] 2 , (7) 



dYd\n(Q 2 /A 2 ) tt y R 2 Q 

where, as before, Y = ln(l/x). The first term on the right-hand-side is the 
growth of g due to DGLAP evolution, and the second is the decrease due to 
gluon-gluon recombination. To gain insight into the origin of the shadowing 
term, note that the number, n, of gluons per unit rapidity interval is xg(x, Q 2 ). 
Moreover the gluon-gluon cross section a(gg) ~ -Ka 2 s /Q 2 and so 

prob. of recomb. ~ - [ x 9? , (8) 
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where nR 2 is the transverse area populated by gluons. The GLR equation 
effectively resums the 'fan' diagrams generated by the branching of QCD 
Pomerons. Recently there has been much activity in this area EJc3, which has 
resulted in an improved knowledge of the triple-Pomeron vertex. 

The structure of the triple-Pomeron vertex can be extracted from an 
equation lB for a quantity, N(r,b,Y), closely related to the cross section for 
the interaction of a qq dipole of transverse size r with the proton target 

a(r,Y) = 2 / d 2 bN{r,b,Y). (9) 

b is the impact parameter of the interaction. In the short-distance approxi- 
mation (r -C b) the non-linear shadowing equation takes the simplified form 

OY 7T I J 

in the large Nq limit, where K is the BFKL kernel and 

f°° d£ 

N(r,b,Y) = J jJ (£r)N(£,b,Y). (11) 

There have been faazp recent attempts to use ([lO]) to calculate the effect 
of gluon shadowing E30. Both give similar results. Fig. 2 shows the results 
of Kimber et al. We see that shadowing is small in the domain accessible to 
experiments at HERA. However almost a factor of two suppression is antic- 
ipated for x ~ few xl0~ 6 . This domain may be accessible at the LHC for 
either prompt photon or high-Qr Drell-Yan production at large rapidities, 
both of which can proceed via the subprocess gq — > jq. Another approach E3 
evaluates the corrections to dF^ jd In Q 2 due to twist-4 gluon recombination. 

If Fig. 2 shows that shadowing is small at HERA, does this rule out the 
saturation models, such as Eil? Not necessarily; it is hard to distinguish the 
Q 2 dependence with shadowing present from the pure DGLAP Q 2 depen- 
dence, since the difference can be removed by an adjustment of the starting 
distributions E3. However a distinction may be possible at the LHC. A much 
better way to identify shadowing is to study the A dependence of scattering 
on nuclei. 

6 Unintegrated parton distributions 

The natural framework with which to discuss DIS and related hard scattering 
processes is to use parton distributions f a (x,k 2 ,n 2 ), unintegrated over the 
transverse momentum kt of the parton, together with the kt factorization 
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Figure 2. The effect of shadowing p« the integrated gluon distribution xg(x,Q 2 ) at Q 2 = 
4 GeV 2 , taken from Kimber et al. C3. 



theorem EL The unintegrated distributions have the advantage that they 
exactly correspond to the quantities which enter the Feynman diagrams and 
therefore allow for the true kinematics of the process even at LO. Here, for 
simplicity, we will discuss the gluon and so omit the subscript a — g. The 
distribution f(x,kf,n 2 ) depends on two hard scales — kt and the scale [i of 
the probe. The scale \i plays a dual role. On the one hand it acts as the 
factorization scale, while on the other hand it controls the angular ordering 
of the gluons emitted in the evolution. The two-scale distribution satisfies the 
CCFM equation Ell, which embodies both DGLAP and BFKL evolution. In 
practice, it is complicated to solve the equation, and, oap to now, it has only 
proved to be practical within Monte Carlo generators E3. 

To gain insight, recall that both DGLAP and BFKL evolution are essen- 
tially equivalent to ordered evolution in the angles of the emitted gluons. In 
the DGLAP collinear approximation the angle increases due to the growth 
of kt, while in BFKL the angle (0 ~ k t /ke) grows due to the decreasing lon- 
gitudinal momentum fraction as we proceed along the emission chain from 
the proton. The factorization scale fi separates the gluons associated with 
emission from different parts of the process, that is from the beam and tar- 
get protons (in pp collisions) and from the hard subprocess. For example, /j, 
separates emissions from the beam (with polar angle 8 < 90°) from those of 
the target (with 9 > 90°), and from the intermediate gluons associated with 
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the hard subprocess. This separation was proved by CCFM and originates 
from the destructive interference of the different emission amplitudes in the 
angular boundary regions. Since the evolution process is essentially controlled 
by one quantity, the emission angle, we should expect to be able to obtain 
the unintegrated gluon distribution /(cc, A;^, yit 2 ) from a single-scale evolution 
equation. Indeed it is possible to accomplish this and to follow an analytic 
approach where the physical assumptions are more evident and where, in prin- 
ciple, the NLO corrections can be included. The key obsepation is that the 
fi dependence enters only at the last step of the evolution E3. 

To illustrate the last-step procedure, we start from the simplified case of 
pure DGLAP evolution for G = xg, where g is the conventional (integrated) 
distribution 



i-A .„ . 

p 

QQ 



(z)G(^,kfjdz - G(x,k 2 ) j zP gg {z)dz 



dG(x,k 2 ) a s 
d\nk 2 ~ 27r ,- ■■ 

(12) 

Suppose that we were to omit the virtual contribution, then the unintegrated 
density would be 

*«=w=ir p " (!)<! (^ , ) t <i3) 

The virtual contributions do not change the kt of the gluon and may be 
resummed to give the survival probability T that the gluon remains untouched 
in the evolution up to the factorization scale /i 

T(fc t , M ) = exp ^- jf ^ H £ A zP gg {z)dz\ , (14) 

as in the Sudakov form factor. Thus the probability to find a gluon with 
transverse momentum kt (which initiates a hard subprocess with factorization 
scale /i) is 

rl-A fx \ 

p gg( z ) 9 {-> k tJ dz 



f(x,kt,fi 2 ) = T(kt,fi) 



as 
2vr 



(15) 



It is at this last step that the unintegrated distribution becomes dependent 
on ji. Angular ordering requires that the cut-off A = kt/(fJ> + kt) in (|l5|), and 
analogously in ( |l4| ) . The last step causes the distribution to increasingly spill 
over into the domain k t > [i as x decreases, as shown by the dotted curves in 
Fig. 3. 

To include lnll/x) effects, (|l^) is replaced by the,_Linified equation of 
Kwiecinski et al. liJ, and a similar procedure followed tJ. The continuous 
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Figure 3. The continuous curves show the k 2 dependence of the unintegrated gluon 
f(x, k 2 , fj, 2 ) for x = 0.1 and 0.0001 at fj, = 10 GeV obtained from a unified BFKL/DGLAP 
study. The d«*ked curves are obtained from DGLAP evolution. The figure is taken from 
Kimber et al. EJ. 



curves in Fig. 3 are obtained, which are not very different from the previ- 
ous DGLAP results. We conclude that the main physical effects come from 
angular ordering in the last step of the evolution, and not from ln(l/x) terms. 

By imposing angular ordering in both the BFKL and DGLAP terms the 
integral up to fi 2 of / does not equal the integrated gluon. An evaluation of 
/ which imposes the_equality, but does not have complete angular ordering, 
has also been made L3. The difference is a NLO effect c3. 



7 Generalized parton distributions 



Two relevant examples, illustrating the need for generalized (or skewed) par- 
ton distributions, are shown in Fig. 4. The first shows either deeply virtual 
Compton scattering (DVCS) or the electroproduction of vector mesons at 
t = 0. The second is double-diffractive Higgs production in pp collisions. 



These processes involve off-diagonal proton matrix elements {p'\ . 
longitudinal components of momentum transfer, x' ^ x. We have 



dcr, » 



7P or Vp) 



t=o 



= dz 



\p), with 



(16) 
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Figure 4. Examples of processes which depend on the generalized (skewed) gluon distribu- 
tion: (a) DVCS and the electroproduction of vector mesons, (b) double-diffractive Higgs 
production. 



where /j, 2 
a(pp 



z(l-z)Q 2 
p + H + p) = 



k 2 



and 



M J 



M 2 H 



f x 2 ,x 2 ;£ 



Ml 



At small x and x' < < f-^P, 
butions are given by Ej 



(17) 

applicable to these processes, the generalized distri- 



f(x,x' ; el[i 2 ) = Rf(xJliJ?) 



R 



d\n£ 2 



(yT(tt,li)xg(x,%)), (18) 



where the survival probability T is given by (|14[)_Fo 
the skewed to the diagonal distribution is known Ea. It leads, for example, to 
an enhancement of R 2 ~ 2 for a(jp — > Tp) at HERA Ell which seems to be 
required by the data, and an enhancement of R A ~ 2 of <r(pp — > p + H + p) 
at the LHC B. 

Data for -f*P ~ > PlP now exist from high energies (HI, ZEUS), through 
intermediate energies (E665,NMC), down to low energies (HERMES). At 
high energies (x < 0.01) the gluon mechanism shown in Fig. 4(a) dominates, 
whereas at low energies Cx. > 0.1) quark exchange takes over. The calculation 
of Vanderhaeghen et al. E3 gives a satisfactory description. 

Good data are now starting to accumulate on the classic DVCS process, 
which allows generalized parton distributions to be studied without the com- 
plications associated with the vector meson wave function. Allowance must 
be made for the Beth£-Heitler process, but the first indications are that the 
theory expectations cB are in good agreement with the data. 

The original aim for observing DVCS was to measure the generalized par- 
ton distributions, and then to determine the angular momentum distributions 
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of the partons in the proton ej. However generalized distributions are inter- 
esting in their own right and it is an active area for theoretical study. There 
are many interesting observable asymmetries. The evolution of the distribu- 
tions is known to NLO. In analogy to the Mellin moments of ordinary partons, 
we should consider the Gegenbauer moments of the generalized distributions. 
The inverse transform is known cil. There hasie been attempts to calculate the 
distributions in the non-perturbative region e3. An area of present activity-is 
the preservation oLgauge invariance, which requires a twist-3 contribution ej. 
The talk by Diehl e3 gives an overview of some of these developments. 



Acknowledgements 

I thank Markus Diehl, Martin Kimber, Jan Kwiecinski, Misha Ryskin and 
Robert Thorne for valuable help with this talk, and Rosario Nania for arrang- 
ing such a splendid Workshop in Bologna. 



References 



1. HI Collaboration : Eur. Phys. J. C13 (2000) 609; |hep-ex/0012052 



10 
11 
12 
13 



1G. 
17. 
18. 



hep-ex/0012053| . 

2. ZEUS Collaboration: |hep-ex/0105090 . 

3. DO Collaboration: Phys. Rev. Lett. 86 (2001) 1707. 

4. CDF Collaboration: Phys. Rev. D64 (2001) 032001. 

5. CCFR/NuTeV Collaboration: Ph ys. Rev. Lett. 
tLep-ex/001000l| |hep-ex/0102049|. 

6. E866 Collaboration: ^iep-ex/0103030 . 



86 (2001) 2742; 



7. M. Botje, Eur. P hys. J. CIO (2000 ) 285 



J. Pumplin et ah, iep-ph/0101032 , see als o D. Stump, these proceedings. 
W.T. Giele, S.A. Keller, D.A. Kosower, |iep-ph/0104052j ; see also S. 
Alekin, Eur. Phys. J. CIO (1999) 395. 
MRST(1999), E ur. Phys. J. C14 (2000) 133. 

D. Stump et al., tiep-ph/0101051, see also D. Stump, these proceedings. 



MRST(2001), see R. Thorne et al., these proceedings ( |hep-ph/0106075D . 
D. Haidt, these proceedings. 

14. K.J. Golec-Biernat, M. Wiisthoff, Phys. Rev. D59 ( 1998) 014017. 

15. V. Barone, C. Pascaud, F. Zomer, [hep-ph/0004268. 



CTEQ Collaboration: Eur. Phys. J. C12 (2000) 375. 
J. Kwiecinski, A.D. Martin, A. Stasto, Phys. Rev. D56 (1997) 3991. 
V.S. Fadin, L.N. Lipatov, Phys. Lett. B429 (1998) 127; G. Camici, M. 
Ciafaloni, Phys. Lett. B430 (1998) 349. 



DCPT-01-58: submitted to World Scientific on February 1, 2008 11 



19. M. Ciafaloni, D. Colferai, G. Salam, Phys. Rev. D60 (1999) 114036 and 
refs. therein; see also G.P. Sala m, Acta. Phys. Polon. B30 (1999) 3679. 

20. R.S. Thorne, [hep-pri/0103210| . 

21. J. Kwiecinski, A.D. Martin, P.J. Sutton, Z. Phys. C71 (1997) 585. 

22. G. Altarelli, R.D. Ball, S. Forte, Nucl. Phys. B599 (2001) 383. 

23. L.V. Gribov, E.M. Levin, M.G. Ryskin, Phys. Rep. 100 (1983) 1. 

24. Ia. Balitsky, Nucl. Phys. B463 (1996) 99; J. Bartels, M. Wusthoff, 
Z. Phys. C66 (1995) 157; A. Mueller, G. Salam, Nucl. Phys. B475 
(1996) 293; R. Venugopalan, Acta. Phys. Polon. B30 (1999) 3731; A. 
Bialas, H. Navel et, R. Peschanski , Nucl. Phys. B593 (2001) 438; E. 
Levin, U. Maor, |hep-ph/0009217| ; H. Weigert et al Phys. Rev. D62 
(2000) 114005; E. Levin, K. Tuchin, [hep-pn/0012167| ; Y. Kovchegov, L. 
McLerran, Phys. Rev. D60 (1999) 054025; and references therein. 

25. Y.V. Kovchegov, Phys. Rev. D60 (1999) 034008; D61 (2000) 074018. 

26. MA. Kimber, J. K wiecinski, A.D. Ma rtin, Phys. Lett. B508 (2001) 58. 

27. M. Lublinsky et al., friep-ph/0 10232 1| . 

28. J. Blumlein et al., Phys. Lett. B504 (2001) 235. 

29. A. Capella et al., Nucl. Phys. B593 (2001) 336. 

30. S. Catani, M. Ciafaloni, F. Hautmann, Nucl. Phys. B366 (1991) 135; 
J.C. Collins, R.K. Ellis, Nucl. Phys. B360 (1991) 3; E.M. Levin et al, 
Sov. J. Nucl. Phys. 54 (1991) 867. 

31. M. Ciafaloni, Nucl. Phys. B296 (1988) 49; S. Catani, F. Fiorani, G. 
Marchesini, Nucl. Phys. B336 (1990) 18. 

32. H. Jung, G.P. Salam, Eur. Phys. J. C19 (2001) 351; H. Kharraziha, L. 
Lonnblad, JHEP 98003 (1998) 006. 

33. M. Kimber, A.D. Martin, M.G. Ryskin, Phys. Rev. D63 (2001), 114027. 

34. M. Kimber et al., Phys. Rev. D62 (2000) 094006. 

35. V.A. Khoze, A.D. Martin, M.G. Ryskin, Eur. Phys. J. C14 (2000) 525. 

36. A.G. Shuvaev et al., Phys. Rev. D60 (1999) 014015. 

37. A.D. Martin, M.G. Ryskin, T. Teubner, Phys. Lett. B454 (1999) 339; 
L. Frankfurt, M. McDermott, M. Strikman, JHEP 9902 (1999) 002. 

38. M. Vanderhaeghen et al., Phys. Rev. Lett. 80 (1998) 5064. 

39. L. Frankfurt, A. Freund, M. Strikman, Phys. Rev. D58 (1998) 114001; 
D59 (1999) 119901E. 

40. X. Ji, Phys. Rev. D55 (1997) 7114. 

41. A.G. Shuvaev, Phys. Rev. D60 (1999) 116005. 

42. V. Petrov et al., Phys. Rev. D57 (1998) 4325. 

43. A. Radyushkin et al., Phys. Rev. D63 (2001) 114012, and refs. therein. 

44. M. Diehl, these proceedings. 



DCPT-01-58: submitted to World Scientific on February 1, 2008 12 



